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ABSTRACT 

We  present  a  geometrical  study  of  a  missile  seeker  illuminating  a  ship  target. 
Modern  seekers  employ  various  methods  to  recognise  their  target  amid  back¬ 
ground  vessels  and  electronic  countermeasures  such  as  chaff  and  active  decoys. 
We  focus  on  one  such  possible  method:  the  seeker’s  measurement  of  a  ship’s 
down-range  extent.  We  aim  to  find  a  representative  value  for  this  extent  given 
that  the  seeker  can  approach  the  ship,  centre  on,  from  any  direction.  This 
representative  value  can  then  be  used  in  studies  of  electronic  attack  in  the  gen¬ 
eration  and  presenting  to  the  seeker  of  a  realistic  false  ship  target.  Ensuring 
the  false  target  has  a  statistically  favourable  range  extent  might  make  it  ap¬ 
pear  more  realistic  to  the  seeker.  On  the  other  hand,  a  lack  of  any  statistically 
dominant  range  extent  implies  that  we  needn’t  select  any  particular  value  to 
present  to  the  seeker.  Either  way,  this  report’s  analysis  gives  guidance  in  what 
range  extent  to  specify  when  constructing  a  false  target. 

The  simple  scenarios  we  consider  turn  out  to  present  no  statistically  dom¬ 
inant  range  extents.  Extrapolating  from  these,  we  conclude  that  it’s  probably 
not  necessary  to  focus  strongly  on  any  one  range  extent  when  constructing  a 
false  target. 
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Statistics  of  Radial  Ship  Extent  as  Seen  by  a  Seeker 

Executive  Summary 

This  report  studies  the  geometry  and  associated  statistics  of  a  missile  seeker  illuminating 
a  ship  target.  The  seeker  must  isolate  its  target  from  other  nearby  sea- going  vessels  that 
it  has  no  interest  in  targeting.  Modern  seekers  employ  various  methods  to  recognise  their 
target  amid  background  vessels  and  electronic  countermeasures  such  as  chaff  and  active 
decoys.  In  this  report  we  focus  on  one  of  these  possible  methods:  the  seeker’s  measurement 
of  a  ship’s  down-range  extent.  A  wide-beam  radar  on  a  seeker  typically  returns  little  or  no 
information  about  a  ship’s  cross-range  extent,  so  down-range  extent  becomes  the  dominant 
extent  parameter  that  it  analyses. 

The  seeker  measures  the  target’s  range  extent  in  some  way  that  we  don’t  need  to  spec¬ 
ify;  we  stipulate  only  that  the  seeker  illuminates  the  target  with  a  beam.  This  beam  could 
be  wide  (from  a  radar),  or  narrow  (from  a  laser).  Our  study  aims  to  find  a  representative 
value  for  the  measured  range  extent,  given  that  the  seeker  can  approach  the  ship,  centre 
on,  from  any  direction. 

To  implement  a  false-target  electronic  attack,  we  require  a  representative  value  of 
range  extent  in  order  to  generate  and  present  to  the  seeker  a  realistic  ship  image.  If 
certain  seeker-measured  range  extents  turn  out  to  be  statistically  favoured  when  all  angles 
of  the  seeker’s  incoming  direction  have  been  examined,  then  creating  our  false  target  with 
one  of  these  values  might  make  it  appear  more  realistic  to  the  seeker.  On  the  other  hand,  if 
no  range  extents  are  statistically  dominant,  we  are  not  forced  into  selecting  any  particular 
value  to  apply  to  the  range  extent  of  the  false  target  that  we  generate.  Either  way,  the 
outcome  of  the  following  analysis  gives  guidance  in  what  range  extent  to  specify  when 
constructing  a  false  target. 

We  find  that  a  wide  seeker  beam  yields  down-range  extent  values  near  to  the  ship’s 
length,  while  a  narrower  beam  yields  extent  values  near  to  the  ship’s  width.  But  no  very 
strongly  dominant  range  extents  emerge  from  our  statistical  analysis.  We  can  expect  the 
seeker  to  be  aware  of  these  statistics  when  it  makes  assumptions  about  the  nature  of  its 
target.  Our  conclusion  is  that  it’s  probably  unnecessary  to  focus  on  any  one  particular 
range  extent  when  constructing  a  false  target  presentable  to  a  seeker. 

The  results  of  this  report  allow  a  selection  of  more  representative  range  extent  values 
for  ongoing  work  into  false-target  generation  for  countering  incoming  missiles. 
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1  Introduction 

This  report  presents  a  study  of  the  geometry  and  associated  statistics  of  a  missile  seeker 
illuminating  a  ship  target.  The  seeker  must  isolate  the  ship  from  a  background  that  it 
has  no  interest  in  targeting,  such  as  other  nearby  sea-going  vessels  and  electronic  counter¬ 
measures  such  as  chaff  and  active  decoys.  In  this  report  we  focus  on  one  of  these  possible 
methods:  the  seeker’s  measurement  of  a  ship’s  down-range  extent.  A  wide-beam  radar 
on  a  seeker  typically  returns  little  or  no  information  about  a  ship’s  cross-range  extent,  so 
down-range  extent  becomes  the  dominant  extent  parameter  that  it  analyses. 

The  seeker  measures  the  target’s  range  extent  in  some  way  that  we  don’t  need  to  spec¬ 
ify;  we  stipulate  only  that  the  seeker  illuminates  the  target  with  a  beam.  This  beam  could 
be  wide  (from  a  radar),  or  narrow  (from  a  laser).  Our  study  aims  to  find  a  representative 
value  for  the  measured  range  extent,  given  that  the  seeker  can  approach  the  ship,  centre 
on,  from  any  direction. 

To  implement  a  false-target  electronic  attack,  we  require  a  representative  value  of  range 
extent  in  order  to  generate  and  present  to  the  seeker  a  realistic  image  of  the  ship  it  seeks.  If 
certain  seeker-measured  range  extents  turn  out  to  be  statistically  favoured  when  all  angles 
of  the  seeker’s  incoming  direction  have  been  examined,  then  creating  our  false  target  with 
one  of  these  values  might  make  it  appear  more  realistic  to  the  seeker.  On  the  other  hand,  if 
no  range  extents  are  statistically  dominant,  we  are  not  forced  into  selecting  any  particular 
value  to  apply  to  the  range  extent  of  the  false  target  that  we  generate.  Either  way,  the 
outcome  of  the  following  analysis  gives  guidance  in  what  range  extent  to  specify  when 
constructing  a  false  target. 

As  an  example  of  relevant  statistics,  a  seeker  will  presumably  not  reject  a  ship-sized 
target,  but  it  might  well  reject  targets  that  are  too  short  in  range  extent,  because  these 
could  be  smaller  vessels  such  as  fishing  boats.  So  a  typical  question  to  address  is  how  likely 
it  is  that  the  range  extent  of  a  real  ship  will  be  measured  to  have  its  smallest  possible  value. 

The  geometrical  analysis  in  this  report  is  aimed  at  complementing  other  analyses  that 
focus  on  the  waveform  aspects  of  false-target  generation.  Most  work  in  this  field  seems  to 
focus  on  the  high-quality  images  produced  by  ISAR.  Eor  example.  Pace  et  al.  [1]  describe  a 
hardware  implementation  of  an  image  synthesiser  for  countering  ISAR.  Despite  their  using 
a  much  higher  granularity  than  our  wide-beam  analysis,  along  with  Doppler  processing, 
still  only  down-range  information  is  considered.  Benedek  [2]  proposes  analysing  the  images 
of  range-Doppler  maps  by  correlating  their  successive  frames;  again,  down-range  informa¬ 
tion  is  being  heavily  relied  upon.  Kostis  [3]  describes  the  production  of  false-target  images 
presentable  to  an  ISAR,  and  focuses  on  creating  as  much  realism  as  possible  through  a 
fine-grained  analysis  of  how  the  ISAR  analyses  target  characteristics. 


2  Ship  Models  Used  and  Assumptions  Made 

We  use  two  ship  models  and  two  beam  models  to  construct  four  scenarios  of  a  ship  be¬ 
ing  illuminated  by  the  radar  beam  of  an  incoming  seeker.  All  ships  and  beams  are  two 
dimensional.  Eirst  we  consider  a  wide  radar  beam,  the  domain  of  current  seekers,  and 
analyse  such  a  beam  illuminating  an  elliptical  ship  and  a  rectangular  ship.  After  this  we 
consider  a  narrow  beam,  which  could  be  the  domain  of  future  seekers,  and  again  perform 
the  analysis  with  this  beam  illuminating  elliptical  and  rectangular  ships. 
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The  ship  models  have  length  L  and  width  IT,  and  without  loss  of  generality  we  always 
assume  W  ^  L.  (Note  that  the  ship’s  width  is  often  called  its  beam,  but  to  avoid  confusion 
with  the  seeker’s  radar  beam,  we  don’t  use  this  term.)  In  all  cases  we  describe  the  beam 
of  the  seeker  as  intersecting  the  centre  of  the  ellipse  or  rectangle  at  some  angle  6  to  the 
“length”  line  running  fore  to  aft  of  the  ship,  measured  clockwise  from  the  forward  direction. 
The  seeker  can  approach  from  any  bearing  with  equal  likelihood,  making  all  values  of  6 
equally  likely. 

At  some  given  moment  the  seeker  measures  the  ship’s  range  extent  to  be  some  length  E. 
Given  that  the  seeker  can  approach  from  any  bearing,  we  wish  to  calculate  the  extent  E 
and  to  examine  its  statistics  as  functions  of  L,W,  and  6: 

-  its  probability  density  p{E), 

-  its  mean  (E),  and 

-  the  cumulative  probability  C{E)  of  measuring  any  extent  up  to  some  given  E. 

The  fore-aft  and  left-right  symmetry  of  the  two  ship  models  ensures  we  need  only 
consider  incoming  seeker  angles  in  the  first  quadrant:  0^0^  90°.  This  is  due  to  a 
combination  of  all  quadrants  being  equally  likely  to  contain  the  seeker,  and  the  fact  that 
the  value  of  E  for  any  value  of  9  is  always  equal  to  the  E  for  the  corresponding  angle 
determined  by  symmetry  in  the  quadrant  0  — >•  90°.  In  other  words,  we  can  replace  any 
value  of  0  >  90°  by  a  corresponding  0  in  the  interval  0  =  0—)-  90°. 

To  construct  (E)  and  C{E)  requires  knowledge  of  p{E).  We  will  begin  with  the 
only  probability  density  available,  p{9)  the  probability  density  for  9,  and  use  it  to  con¬ 
struct  p{E).  (It’s  conventionally  understood  that  these  are  two  different  functions  “p”.)  By 
definition  of  the  probability  density,  the  probability  that  a  seeker  will  be  inbound  in  the 
angle  interval  0  — )•  0  -|-  d0  is  p{9)  |d0|.  As  discussed  above,  we  will  treat  the  first  quadrant 
as  being  the  only  one  that  the  seeker  accesses,  in  which  case  its  probability  of  being  found 
there  is  one: 

/2 

p{9)d9  =  l.  (2.1) 

But  p{9)  is  independent  of  9,  so  we  conclude  that  p{9)  =  2/7r. 

If  the  angle  interval  9  ^  9  +  d9  were  to  map  uniquely  to  a  corresponding  extent  in¬ 
terval  E  ^  E  +  dE,  the  extent  probability  in  the  extent  interval  would  equal  the  angle 
probability  in  the  angle  interval: 


p{E)\dE\=2/7r\d9\,  (2.2) 

since  clearly  the  fraction  of  measurements  p{E)  IdEl  falling  into  an  extent  bin  must  be  the 
same  as  the  fraction  2/7r  |d0|  falling  into  the  corresponding  angle  bin — because  the  same 
measurements  occupy  both  bins;  the  two  bins  are  just  different  ways  of  labelling  the  same 
measurements.  We’ll  find  this  to  be  the  case  for  both  ellipse  models,  but  for  the  rectangle 
models  there  can  be  two  angle  bins  that  contribute  to  one  extent  bin.  We’ll  encounter 
those  in  due  course  and  will  modify  (2.2)  to  sum  the  probabilities  in  two  angle  bins  instead 
of  one. 

In  the  next  sections  we  determine  the  extent  E  and  its  statistics  p{E),  (E),  and  C{E) 
for  each  of  the  four  scenarios  described  on  the  preceding  page. 
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Figure  1:  Definitions  of  terms  used  in  this  analysis 

3  Scenario  1:  Ellipse  Illuminated  by  a 

Wide  Beam 

We  model  the  ship  as  an  ellipse  in  two  dimensions,  defining  its  extent  E  in  Figure  1. 

3.1  The  Extent  E 

To  calculate  the  extent  E  as  a  function  of  0,  begin  by  noting  that  the  slope  of  the  seeker’s 
line  of  sight  of  the  ship  at  angle  9  relative  to  the  fore-aft  axis  (the  y  axis  in  Figure  1) 
is  cot  9]  hence  the  slope  of  the  rightmost  dotted  perpendicular  line  in  that  figure  is  —tan  9. 
This  dotted  line  therefore  has  equation  y  =  yo  —  tan  9  {x  —  x^).  The  requirement  that  this 
line  intersect  the  ellipse  at  a  single  point  determines  (a:o,yo)i  the  extent  E  is  then 

given  by  (E/2)^  =  Xq  +  y^.  The  ellipse  has  equation  ~  which  becomes 

y  =  Ly^l/d  —  x"^  '  in  the  first  quadrant.  So  we  require  that  the  two  equations,  of  the 

dotted  line  through  (xo,2/o)  of  the  ellipse,  have  one  solution  {x,y): 

y  =  Vo  —  tan  9  {x  —  Xq)  (dotted  line) , 
y  =  (ellipse)  .  (3.1) 

Equating  the  right-hand  sides  of  (3.1)  gives  (with  t  =  tan0  =  Xq/j/o  for  shorthand) 

Vo  -  t{x  -  Vot)  =  -  .  (3.2) 

Setting  k  =  1  + 1'^  =  seE  9  and  squaring  both  sides  of  (3.2)  gives  a  quadratic  in  x: 

(t^  +  x^  -  2yQktx  -  ^  =  0 ,  (3.3) 
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Figure  2:  Left:  The  extent  in  units  of  width  for  a  length-to-width  ra¬ 
tio  L/W  =  5,  where  it’s  sufficient  to  plot  only  angles  in  the  first  quadrant. 
Right:  A  polar  plot  of  the  same  extent  in  all  four  quadrants.  The  vertical 
axis  here  runs  along  the  length  of  the  ship. 


which  will  have  a  unique  solution  (corresponding  to  the  dotted  line  intersecting  the  ellipse 
at  a  single  point)  if  its  discriminant  “6^  —  4ac”  equals  zero;  that  is,  if 


4y2feV-4(t^  +  ^)(yo"fc^-:-)  =0. 


L^\ 


W^J 


.2;  2 


L^\ 


4  / 


Equation  (3.4)  simplihes  to 


2  ^ 


But  the  extent  is  given  by  {E(2)'^  =  Xq  +  j/q,  or 

=  4(xo  +  Vo)  =  4(t^yo  +  Vo)  =  4A:yo  • 

Now  substitute  the  expression  for  yg  from  (3.5)  into  the  last  term  of  (3.6),  giving 


^  “T- 


+ 


sec^  0 


(3.4) 


(3.5) 


(3.6) 


(3.7) 


which  simplifies  to 

E  =  ^/ sin^  9  +  L'^  cos^  9 ' . 


(3.8) 


The  extent  in  units  of  width,  E/W,  is  plotted  as  a  function  of  6  in  Figure  2  for  a  length- 
to-width  ratio  L/W  =  5. 


3.2  Probability  Density  p{E) 


Equation  (2.2)  requires  that  we  relate  d0  to  dE.  Differentiating  (3.8)  with  respect  to  9 
leads  to 


^  _ _ E _ 

dE  {W"^  —  L'^)  sin  9  cos  9 


(3.9) 
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Figure  3:  Probability  density  p{E)  from  (3.13)  for  a  length-to-width  ratio 
L/W  =  5 


This  expression  is  always  negative  in  the  interval  0  =  0  — )■  7r/2,  meaning  E  decreases  while 
6  increases,  which  is  expected  and  evident  in  Figure  2.  This  one-to-one  mapping  of  angle 
to  extent  implies  that  every  angle  bin  0  ^  9  +  dO  corresponds  uniquely  to  one  extent 
bin  E  ^  E  +  dE,  as  in  (2.2).  We  conclude  that  (2.2)  can  be  applied  as 


TT  dE 
2E 

—  W^)  sin  9  cos  9 


p{E)  =  - 

TT 


d9 


dE 


(3.10) 


To  express  p{E)  wholly  in  terms  of  E  we  need  to  write  sin  9  and  cos  9  in  terms  of  E.  Begin 
by  writing  (3.8)  as 

E"^  =  W'^W9  +  L^{l-W9),  (3.11) 


which  leads  to 


sin0 


L‘^-E^'  ^  lE^  -  ' 

L^-W^  ’  ^  -  V 


(3.12) 


in  the  first  quadrant.  These  convert  (3.10)  to  an  expression  wholly  in  terms  of  E: 


p{E) 


2E 

7rv^(L2-E2)(E2-iy2)'' 


(3.13) 


The  probability  density  p{E)  for  a  length-to-width  ratio  5  is  plotted  in  Figure  3.  Although 
(3.13)  shows  that  p{E)  goes  to  infinity  when  the  extent  equals  either  the  ship’s  width  or 
its  length,  the  divergence  at  these  extremes  is  not  strong  enough  to  prevent  p{E)  from 
integrating  to  1  between  these  two  values,  as  is  required  of  a  probability  density.  The 
high  density  near  the  extreme  values  of  W  and  L  is  due  to  the  change  in  extent  being 
of  a  smaller  order  as  the  angle  9  sweeps  through  an  infinitesimal  range  near  these  values, 
as  compared  with  the  change  in  extent  for  other  angle  values.  It  can  be  likened  to  the 
comparatively  long  time  that  a  piston  spends  at  the  top  and  bottom  of  its  travel. 

This  bi-modality  of  p{E)  is  its  key  feature,  and  it  implies  that  p{E)  cannot  simply 
be  rendered  down  to  a  single  representative  value.  (Contrast  this  with  a  bell-shaped 
probability  density,  which  can  be  characterised  by  its  mean,  together  with  a  spread  about 
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Figure  4-  The  mean  extent  (E)  in  units  of  width  W  for  a  range  of  length- 
to-width  ratios  a,  found  by  evaluating  (3.15)  numerically 


that  mean.)  So  while  we  might  choose  to  represent  p{E)  by  the  resulting  mean  extent,  (E), 
we  must  remember  that  this  value  is  representative  only  in  a  limited  context. 


3.3  Mean  Extent  {E) 

Refer  to  (3.13)  to  write 

- C ^  f C  ■ 

Write  this  integral  dimensionlessly  by  setting  x  =  EjW  and  a  =  LfW ,  giving 


{E)  _  2  r  x^dx 

tE  TT  J I  ^ {p?-  —  X^)(X^  —  1)  ' 


The  result  of  the  integration  is  an  elliptic  integral  of  the  second  kind: 

(E)  2a  f  TT  —  1  \ 

2  ’  j  ’ 


(3.16) 


where  the  “E”  on  the  right-hand  side  of  (3.16)  is  the  elliptic  integral  function,  not  the 
extent[4].  In  fact,  the  curve  of  (E) /W  versus  a  is  more  or  less  straight  for  all  reasonable 
values  of  a,  as  shown  in  Figure  4.  This  analysis  is  found  to  agree  well  with  a  Monte  Carlo 
simulation  in  which  (E)  results  from  averaging  values  of  E  that  are  calculated  from  (3.8) 
for  several  thousand  values  of  0  chosen  uniform  randomly. 

Best-fitting  the  curve  in  Figure  4  with  a  line  over  a  range  of  values  of  a  gives 


{E) /W  ~  0.6q;  -I-  0.4 . 


(3.17) 


This  differs  little  from  an  approximation  (E) /W  ~  O.Sa -|- 0.5,  meaning  that  (E)  is  ap¬ 
proximately  just  the  average  of  the  ship’s  length  and  width. 
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Figure  5:  C{E)  for  a  length-to-width  ratio  5,  found  by  evaluating  (3.19) 
numerically  with  a  =  5 


3.4  Cumulative  Probability  of  the  Extent,  C{E) 


The  cumulative  probability  of  the  extent  is  the  probability  that  the  extent  will  be  measured 
to  have  some  value  between  the  minimum  possible,  the  width,  and  a  given  value  E: 


C{E)  =  [""piE) 

Jw 


dE 


(3.13)  2 


EdE 


Jw  V(L2  -E2)(E2  -  WE 

As  in  the  previous  section,  set  x  =  E/W  and  a  =  L/W  to  write 


(3.18) 


(3.19) 


C{E)  has  a  reasonably  fixed  form  for  length-to-width  ratios  of  a  =  5, 10,  and  15.  The  case 
of  a  =  5  is  shown  in  Figure  5. 

As  a  concrete  example  we  calculate  C{{E)),  the  probability  of  measuring  any  extent 
less  than  the  mean  extent,  for  length-to-width  ratios  of  a  =  5, 10,  and  15. 

Length-to-width  ratio  a  =  5:  Equation  (3.15)  with  a  =  5  gives  (E) /W  ~  3.34.  [Com¬ 
pare  this  to  the  simpler  estimate  (3.17)  which  gives  a  value  of  3.4.]  Now  set  E/W  to  3.34 
in  (3.19)  and  evaluate  that  integral  numerically,  giving  C{{E))  ~  0.451. 

As  a  sensitivity  comparison,  set  E/W  to  the  simpler  estimate  of  3.4  and  evaluate  (3.19)  numerically 
to  find  C{(E))  ~  0.462.  Also,  setting  E/W  equal  to  the  mean  of  length  and  width,  3  [according 
to  the  comment  just  after  (3.17)]  and  evaluating  (3.19)  numerically  gives  C{(E))  ~  0.392. 

Length-to-width  ratio  a  =  10:  Equation  (3.15)  with  a  =  10  gives  {E)/W  ~  6.47.  [Com¬ 
pare  this  to  the  simpler  estimate  (3.17)  which  gives  a  value  of  6.4.]  Now  set  E/W  to  6.47 
in  (3.19)  and  evaluate  that  integral  numerically,  giving  C{{E))  ~  0.444. 

Length-to-width  ratio  a  =  15:  Equation  (3.15)  with  a  =  15  gives  {E) /W  ~  9.63.  [Com¬ 
pare  this  to  the  simpler  estimate  (3.17)  which  gives  a  value  of  9.4.]  Now  set  E/W  to  9.63 
in  (3.19)  and  evaluate  that  integral  numerically,  giving  C{{E))  ~  0.442. 
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Figure  6:  Definitions  of  terms  used  in  this  analysis 


These  three  cumulative  probabilities  C{{E))  =  0.451,0.444,0.442  are  all  fairly  similar, 
and  we  conclude  that  for  a  wide  range  of  length-to-width  ratios,  on  average  about  44%  of 
the  extent  measurements  will  yield  values  between  the  minimum  possible  (the  width)  and 
the  mean  extent. 


4  Scenario  2:  Rectangle  Illuminated  by  a 

Wide  Beam 

4.1  The  Extent  E 

We  now  replace  the  ellipse  of  Scenario  1  with  a  rectangle  of  width  W  and  length  L  W, 
again  illuminated  by  a  wide  beam,  shown  in  Figure  6.  For  shorthand,  define  D  as  the 
diagonal  length  of  the  (full)  rectangle: 

D  =  VIU2  +  L2'.  (4.1) 

Two  key  angles  will  be  needed  in  our  analysis:  9^),  the  angle  where  the  extent  attains 
its  maximum  D,  and  0^,  the  nonzero  angle  less  than  90°  where  the  extent  has  value  L 
(shown  in  Figure  7).  Figure  6  shows  that  for  0  ^  9  ^  Od  the  semi-extent  is  given  by 
E/2  =  W/2  sm9  +  L/2  cos9.  We  can  treat  the  case  of  9^  <  9  ^  7r/2  using  Figure  6  by 
swapping  L  and  W  and  setting  9  — )•  tt/2  —  9,  which  yields  the  same  expression  for  E.  So 
we  conclude  that  in  the  whole  first  quadrant, 


E  =  W  sin  9  +  L  cos  9  . 


(4.2) 


This  is  shown  schematically  in  Figure  7  and  quantitatively  in  Figure  8. 
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Figure  7;  A  schematic  of  the  extent  E  as  a  function  of  angle  9  from  (4.2), 
with  angle  contributions  to  the  interval  E  — )•  E  +  dE  shown  in  grey 


Figure  8:  A  polar  plot  of  the  four-quadrant  extent  for  L/W  =  5  from  (4.2); 
the  vertical  axis  here  runs  along  the  length  of  the  ship 
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4.2  Probability  Density  p{E) 

The  calculation  of  p{E)  is  now  more  involved  than  in  the  ellipse  model.  For  6  < 

(E  >  L),  there  are  two  angle-bin  contributions  to  p{E)  IdEl,  because  the  same  value  of  E 
can  result  from  either  of  two  angles.  In  that  case  we  add  the  probabilities  for  the  two 
angle  contributions,  which  is  simply  the  requirement  that  the  number  of  “trials”  p{E)  IdEl 
in  the  extent  bin  be  the  sum  of  the  numbers  of  trials  in  the  two  angle  bins  that  contribute 
to  that  extent  bin.  The  analysis  becomes  more  complicated,  but  this  is  purely  due  to 
the  extra  bookkeeping  required  to  keep  track  of  the  two  angle  contributions  to  the  extent 
probability.  The  two  angle  intervals  that  contribute  to  a  given  extent  interval  are  shown 
as  grey  regions  in  Figure  7. 

The  angle  Ojj  =  tan“^  ^  marks  the  greatest  possible  value  of  the  extent,  being  along 
the  diagonal  of  the  rectangle.  The  angle  marks  a  dividing  line:  to  its  left  in  Figure  7  is 
the  region  where  there  are  two  contributions  to  p{E)  (where  the  function  E{9)  is  not  “one 
to  one”),  and  to  its  right  is  the  region  where  there  is  only  contribution  (where  E{9)  is  one 
to  one).  We  must  calculate  p{E)  in  each  region  separately. 

Note  that  (4.2)  gives 

/W  L 

E  =  W  sm9  +  L  cos  9  =  D  I  —  sm9  +  —  cos  9 

=  i7(sin  sin  9  +  cos  9^)  cos  9) 

=  Dcos{9-9d),  (4.3) 

which  is  a  cosine  shifted  to  the  right  hy  9 d-  Hence  E{9)  is  symmetrical  about  9^),  so  we 
can  see  from  Figure  7  that  9^  =  29^.  Equation  (4.2)  also  produces 

\dE/d9\  =  Vd^-E^'  (4.4) 

over  the  whole  range  of  0,  a  result  that  will  be  required  shortly. 


4.2.1  The  Region  W  ^  E  ^  L 


Figure  7  shows  that  this  extent  interval  corresponds  to  the  angle  interval  9^  ^  9  ^  7r/2. 
There  is  a  single  contribution  to  p{E)  |dE|  here  of  2/7r  |d0|.  It  follows  from  (4.4)  that 


and  therefore 


p{E)  \dE\  =  -  |d0| 
vr 


2  |dE| 


p{E) 


2 

7rVD2  -E2" 


(4.5) 

(4.6) 


4.2.2  The  Region  L  <  E  ^  D 

Figure  7  shows  that  this  extent  interval  corresponds  to  the  angle  interval  0  ^  9  <  9]^. 
The  probability  p{E)  \dE\  now  comes  from  two  angle  contributions  (two  angle  bins)  each 
with  probability  p{9)\d9\  =  2/7r|d0|.  In  both  bins  (4.4)  holds,  so  that  each  gives  the 
contribution  on  the  right-hand  side  of  (4.5);  hence  these  bins  add  to  give 

p{E)  |dE|  =  -  , ,  (4.7) 

'  7rVD2-E2  ’  ^  ^ 


10 


UNCLASSIFIED 


UNCLASSIFIED 


DSTO-TN-1220 


Extent  in  units  of  width,  E/W 

Length-to-width  ratio  3 


Extent  in  units  of  width,  E/W 

Length-to-width  ratio  10 


Figure  9:  Probability  densities  for  two  length-to-width  ratios 


so  that 


p{E) 


4 

7rVD2  -E2'’ 


The  probability  density  for  all  values  of  9  is  then 


(4.8) 


p{E) 


IttUd'-e"'’  L  <  E 


(4.9) 


This  is  plotted  for  two  length-to-width  ratios  in  Figure  9. 


4.3  Mean  Extent  {E),  and  Spread  of  the  Extent  cr| 

The  mean  extent  is 


(E)  =  [  Ep{E)dE 

Jw 

_2  EdE  4  EdE 

Jw  VD2  -E2'  Jl  VD^-E^' 


TT 


-\/D2  -  E2 


H — 

J  W  TT 


-\/d2 


E2 


D 

L 


which  simplifies  to 


(E)  =  2/tt  {W  +  L). 


(4.10) 


(4.11) 


This  is  slightly  more  than  the  mean  of  W  and  L,  and  is  plotted  in  Figure  10. 

As  an  exercise  in  calculating  the  standard  deviation  of  E,  consider  that  the  variance 
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Figure  1 0:  The  mean  extent  (E)  in  units  of  width  W  for  a  range  of  length- 
to-width  ratios  a,  found  from  (4.11) 


is  cj|  =  (E^)  —  (E)^.  We  require: 

Jw  ^  Jw 


rD 


B  + 


dE 


VD2-E2'  TT  Jl  VD2  -E2' 


2 

TT 


E 


D2 


E 


- \/ 1)2  —  E2  '  H - sin  ^  — 

2  2  D 


"  4 

+  - 
£;=iy  ^ 


E  /-5 - ^  D2  .  E' 

2^  2  D 


D2  2IUL 
“1“  +  ^ 


D 

E=L 

(4.12) 


Now  we  can  write 


2  /r^2\  /rn\2  ^lUL  4  3 

,2_/c.2\  /P\2_  , - + 

2  TT  TT^ 


4  =  (e2)  -  (E)2  = 


vr^ 


1U2  +  +  21UL 


TT  —  4 
7r2/2  -4_ 


~0.1(L-1U)2, 


(4.13) 


and  the  standard  deviation  of  the  extent  is  ~  0.3  {L  —  W),  remembering  that  L  ^  W. 
The  extent  is  not  peaked  close  to  its  mean;  for  example,  if  IT  =  1  and  L  =  10,  then 
(E)  ~  7  and  ~  3 — and  of  course  there  is  no  peak  at  all  around  7  in  the  right-hand  plot 
in  Figure  9. 


4.4  Cumulative  Probability  of  the  Extent,  C{E) 

The  cumulative  probability  of  measuring  any  extent  up  to  some  value  E  is 


C{E) 


'w 


p{E)  dE 


(4.9) 


2  fE  dE 

TT  Jw  •yo^—E'^'  ’ 

2  fL  dE  _L  4  dE 

TT  Jw  y'£)2_£;2'  '  IT  Jl  y'D^—E^'  ’ 


E^L 

L  <  E  if  D 
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Extent  in  units  of  width,  E/W 


Figure  11:  C{E)  for  a  length-to-width  ratio  1.5,  from  (4.15)  with  L  =  1.5W 


Extent  in  units  of  width,  E/W 
Length-to-width  ratio  5 


Length-to-width  ratio  10 


Figure  12:  Cumulative  probabilities  for  realistic  length-to-width  ratios, 
from  (4.15).  The  rightmost  part  of  the  curve  becomes  less  significant  for 
higher  ratios 


= 


/ 

2 

.  -| 

E 

TT 

Sin 

D 

< 

2 

.  1 

e' 

TT 

\ 

sin 

D 

1  E 


Jfy 


+  1 


sin-i  § 


1  ® 
-  L 


EifL 
L  <  E  ^  D 


(4.14) 


leading  to 


C{E) 


|(sin-i®-0^),  W^E^L 
i  sin-1  E  _  1  ^  L<E^D.' 


(4.15) 


where  Od  =  tan  ^  ^  as  before.  The  cumulative  probability  is  plotted  in  Figure  11  for  a 
length-to-width  ratio  1.5.  More  realistic  examples  are  given  in  Figure  12. 

As  we  did  for  the  elliptical  ship,  we  again  calculate  C{{E)),  the  cumulative  probability 
of  measuring  an  extent  up  to  the  mean  (E).  Applying  (4.15)  requires  us  to  know  whether 
(E)  is  greater  than  or  less  than  L.  Use  (4.11)  to  write,  after  some  short  algebra. 


{E)<L 


L  > 


W 
-  1 


~  1.75  IF. 


(4.16) 
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(So  (E)  will  be  less  than  L  for  all  practical  length-to-width  ratios,  which  is  reasonable.) 
The  required  cumulative  probability  becomes 


cm)  = 


sm 


-1  (E) 


D 


-e 


D  ! 


■  sm 


-1  {E) 


-1, 


L  > 
L  C 


w 

7r/2-l 

W 


~  1.75  IF 
1.75  IF. 


(4.17) 


D  7r/2-l 

Equation  (4.17)  requires  {E)/D,  which  is  given  by  (4.11)  in  terms  of  a  =  L/W  ^  1  as 


(E)  |(«  +  1) 

D  \/a^  +  1' 


(4.18) 


We  can  use  (4.17)  and  (4.18)  to  calculate  C{{E))  for  representative  values  of  a: 


a  =  5  ^  cm)  -  0.41, 
a  =  10  ^  cm)  -  0.43, 
a  =  15  ^  cm)  -  0.43. 


(4.19) 


These  values  don’t  differ  markedly  from  those  calculated  for  an  elliptical  ship  on  page  7: 
there  the  corresponding  values  were  C{{E))  ~  0.45,  0.44,  0.44. 


5  Scenario  3:  Ellipse  Illuminated  by  a 
Narrow  Beam 

In  this  section  and  the  next,  we  again  use  elliptical  and  rectangular  ship  models  in  two  di¬ 
mensions,  but  with  the  seeker  now  illuminating  them  using  a  narrow  beam.  The  elliptical- 
ship  scenario  is  shown  in  Figure  13. 


5.1  The  Extent  E 


Again  it’s  sufficient  to  deal  with  incoming  missile  angles  in  the  first  quadrant.  The  extent 
is  given  by  E  =  2-\/ Xq  +  ',  where  the  point  {xQ^yo)  =  (yotan0,yo)  lies  on  the  ellipse. 
This  point  must  satisfy  the  equation  for  the  ellipse, 

in  which  case  simple  substitution  leads  to 


_  LIF/2 
~  \/L2tan^0  +  IF2'’ 


Xq  =  J/o  f  ^  • 


(5.2) 


Substituting  these  into  E  =  2^/ x‘^  +  y^'  gives  the  extent: 


This  is  plotted  for  the  case  of  L/W  =  5  in  Figure  14. 


(5.3) 


14 


UNCLASSIFIED 


Extent  in  units  of  width,  E/W 


UNCLASSIFIED 


DSTO-TN-1220 


Figure  13:  Definitions  of  terms  used  in  this  analysis 


Figure  I4:  Left:  The  extent  in  units  of  width  for  a  length-to-width  ra¬ 
tio  L/W  =  5,  where  it’s  sufficient  to  plot  only  angles  in  the  first  quadrant. 
Right:  A  polar  plot  of  the  same  extent  in  all  four  quadrants.  Not  surpris¬ 
ingly  this  is  just  the  outline  of  the  ship,  albeit  twice  as  long  and  wide  as  the 
actual  ship,  due  to  the  way  the  polar  plot  presents  the  data. 
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Figure  15:  Probability  density  p{E)  from  (5.6)  for  a  length-to-width  ratio 
L/W  =  5 


5.2  Probability  Density  p{E) 


Because  E  is  always  decreasing  as  0  increases,  only  one  angle  bin  contributes  to  any  range 
bin.  Hence,  using 


Ee 


^^{L^-E^){E^-WE, 


(5.4) 


we  can  write 


p{E)\dE\  =p{e)\de\ 


2 _ ^1^ _ 

^T3vViE-E^){E^-WE' 


(5.5) 


from  which  it  follows  that 


p{E) 


_ 2LW _ 

TvE^iL"^  -  E2)(E2  -  W‘^Y  ' 


(5.6) 


This  certainly  integrates  numerically  to  one  for  various  random  choices  of  L  and  W .  It’s 
plotted  in  Figure  15  for  the  case  of  L/W  =  5. 


5.3  Mean  Extent  {E) 

The  mean  of  E  is 


(E)  =  [  Ep{E)dE 

Jw 


As  before,  set  x  =  E/W  and  a  = 


(5.6)  2LW  E  dE 

vr  Jw  EiE  -  E^){E^  -  WE  ' 

L/W  to  write  this  integral  dimensionlessly: 


(E)  _  2a  r  dx 

tF  J I  E {E  —  x'^){x‘^  —  1)  ' 


(5.7) 


(5.8) 


The  mean  extent  is  plotted  in  Figure  16  for  practical  values  of  a.  As  can  be  seen,  it  ranges 
typically  from  about  1.5  times  the  width  for  shorter  vessels  to  about  2.7  times  the  width  for 
the  longest  ones.  Why  is  it  comparatively  small  and  insensitive  to  the  ship’s  length?  This 
is  due  to  the  use  of  a  narrow  beam  by  the  seeker.  Such  a  beam  is  only  capable  of  sampling 
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Figure  1 6:  The  mean  extent  (E)  in  units  of  width  W  for  a  range  of  length- 
to-width  ratios  a,  found  by  evaluating  (5.8)  numerically 


(or  illuminating)  approximately  the  full  length  of  the  ship  for  a  small  range  of  incoming 
seeker  angles,  which  makes  these  larger  extent  values  comparatively  rare,  in  contrast  to 
the  wide-beam  case.  This  is  also  shown  by  the  probability  density  in  Figure  15,  which 
gives  the  larger  extents  a  lower  probability  density.  There  is  certainly  a  density  peak  at 
an  extent  equal  to  the  ship’s  length,  but  this  peak  is  too  narrow  to  significantly  modify 
the  extent’s  mean. 


5.4  Cumulative  Probability  of  the  Extent,  C{E) 

The  cumulative  probability  of  measuring  any  extent  up  to  some  value  E  is 


CiE)=  [''piE)dE  ^  [ 

Jw  Jv\ 


2LWdE 


w  ttE^{E  -  E^){E^  -W^y 


and  when  again  we  write  x  =  E/W  and  a  =  L/W,  this  becomes 


(5.9) 


(5.10) 


Figure  17  plots  C{E)  for  a  =  5.  Mean  extents  and  cumulative  probabilities  for  represen¬ 
tative  values  of  a  are 


a  =  5 
a  =  10 
a  =  15 


{E)  ~  1.92 
{E)  ~  2.35 
{E)  ~  2.61 


and  C{{E))  ~  0.67, 
and  C((F))~0.73, 
and  C{{E))  ~  0.75. 


(5.11) 


UNCLASSIFIED 


17 


DSTO-TN-1220 


UNCLASSIFIED 


Figure  17:  C{E)  for  a  length-to-width  ratio  5,  found  by  evaluating  (5.10) 
numerically  with  a  =  5 


6  Scenario  4:  Rectangle  Illuminated  by  a 

Narrow  Beam 

In  this  section  we  model  the  ship  as  a  rectangle  illuminated  with  a  narrow  beam.  The 
scenario  is  shown  in  Figure  18  with  the  plot  of  the  extent  shown  in  Figure  19,  which  we’ll 
derive  shortly. 


6.1  The  Extent  E 

Just  as  in  the  previous  rectangle  model  of  Section  4,  0j^  marks  the  greatest  possible  value 
of  the  extent,  and  points  along  the  diagonal  of  the  rectangle: 

1  W 

Ojj  =  tan  —  .  (6.1) 

1j 

Unlike  the  expressions  for  the  extent  E  for  the  previous  two  ship/illumination  models  in 
(3.8)  and  (4.2),  the  current  extent  must  be  specified  separately  for  angles  less  than  9^ 
(“Zl-region  1”)  and  greater  than  6j^  (“Zl-region  2”) 

The  extent  of  the  narrow  beam  across  the  rectangle  is  easily  found  by  inspecting 
Figure  18: 


1  L/  cos  9  , 

0  ^  ^  (D-region  1) 

E=  { 

\W /  sin0 , 

9^)  ^  9  ^  -K /2  (D-region  2) . 

This  is  plotted  in  Figure  20;  not  surprisingly,  the  polar  plot  simply  gives  a  picture  of  the 
rectangle!  It’s  one  corner  of  this  rectangle  that  is  being  shown  in  Figure  19;  of  course,  the 
cartesian  nature  of  Figure  19  warps  the  corner,  as  compared  to  the  right-angled  version  in 
the  polar  plot  of  Figure  20. 
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Figure  18:  Definitions  of  terms  used  in  this  analysis 


Figure  19:  Schematic  of  the  extent  E  as  a  function  of  angle  6  from  (6.2), 
with  angle  contributions  to  the  interval  E  ^  E  +  dE  shown  in  grey  for 
L-region  2 
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Figure  20:  A  polar  plot  of  the  four-quadrant  extent  (6.2)  for  L/W  =  5;  the 
vertical  axis  here  runs  along  the  length  of  the  ship.  As  in  Figure  If,  this  plot 
shows  the  outline  of  the  ship. 


6.2  Probability  Density  p{E) 

Here  we  will  need 


E  tan  9 
E  cot  9 


EVE^-L^'/L, 

E^E'^  -W^'/W , 


D-region  1 
D-region  2 . 


(6.3) 


Also  as  in  the  previous  rectangle  model,  an  angle  9^  marks  a  dividing  line:  to  its  right 
in  Figure  19  is  “L-region  1”  where  there  is  only  one  angle  contribution  to  p{E),  and  to 
its  left  is  “L-region  2”  where  there  are  two  angle  contributions  to  p{E).  Note  that  the  two 
L-regions  don’t  coincide  with  the  two  D-regions. 


9l  =  sin"^ 


IF 

T’ 


compare  (6.1):  9^  =  tan  ^ 


W 

T  ■ 


(6.4) 


We  must  construct  p{E)  separately  for  each  L-region. 


6.2.1  L-region  1:  W  ^  E  ^  L 

Figure  19  shows  that  in  this  region  there  is  a  single  contribution  of  2/tt  |d0|  to  p{E)  \dE\, 
coming  from  D-region  2  in  (6.3): 

p{E)  \dE\  =  -  |d6l|£,_region  2  =  “  ^^^2  _  ^ 


6.2.2  L-region  2:  L  <  E  ^  D 

Figure  19  shows  that  the  probability  p{E)  \dE\  comes  from  two  contributions  of  2/7r  |d0|, 
one  each  from  D-regions  1  and  2.  The  contribution  from  D-region  1  in  (6.3)  is 


-region  1 


2  LldLl 
n  eVL2  -  L^  ’ 


(6.6) 
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Figure  21:  Probability  density  p{E)  of  ship  extent  from  (6.9)  for  a  ship  of 
length-to-width  ratio  1.5,  clipped  at  the  top  to  display  well 


and  the  contribution  from  D-region  2  in  (6.3)  is 


-region  2 


2  W\dE\ 

TT  Ey/E"^  -  W^'  ' 


Adding  these  two  contributions  (6.6)  and  (6.7)  gives 


p(E)  |dE| 


2 

TT 


L 


+ 


W 


eVe^-l^'  eVe2-w^' 


|dE| 


Equations  (6.5)  and  (6.8)  give  the  complete  probability  density  as 


(6.7) 


(6.8) 


p(E) 


w 

TT  EVE'^-W'^' 


+ 


w 


EVE'-^-L'-^'  ^  EVE'-^-W^' 


lU^  E^L 
L  <  E  ^  D 


(6.9) 


A  plot  of  p(E)  versus  E  is  shown  in  Figure  21  for  a  ship  of  length  1.5  times  its  width.  (This 
unrealistic  ratio  is  chosen  to  better  highlight  the  details  of  the  function.)  The  probability 
density  goes  to  infinity  when  the  extent  equals  either  the  ship’s  width  or  its  length.  The 
high  density  near  these  two  values  is  caused  by  the  change  in  extent  being  of  a  smaller 
order  as  the  angle  0  sweeps  through  an  infinitesimal  range  near  each  of  these  two  values, 
as  compared  with  what  0  sweeps  through  near  other  extent  values.  That  is,  the  extent 
change  is  negligible  compared  to  the  angle  change  for  angles  in  the  vicinity  of  the  ship’s 
width  and  length,  which  has  the  effect  that  the  probability  in  the  extent  bins  around  the 
width  and  length  includes  many  more  angles  and  so  must  be  correspondingly  large.  Note 
that  because  p{E)  is  bi-modal,  the  comment  at  the  end  of  Section  3.2  applies  to  its  mean 
value,  which  we  calculate  shortly. 

More  realistic  probability  densities  are  plotted  in  Figure  22,  using  length-to-width 
ratios  of  5  and  10  respectively. 

Note  that  the  entire  density  still  integrates  to  1  as  required:  equation  (6.24)  ahead — 
the  cumulative  probability  C{E) — returns  a  value  of  1  when  evaluated  at  the  largest  value 
of  extent,  E  =  D. 
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Length-to-width  ratio  5  Length-to-width  ratio  10 

Figure  22:  Probability  densities  from  (6.9)  for  more  realistie  length-to-width 
ratios,  showing  that  the  rightmost  peak  in  the  density  has  less  support  at 
higher  ratios 


6.3  Mean  Extent  (E) 

The  mean  of  E  is 


{E)=  [  Ep{E)dE. 

Jw 


Use  (6.9)  to  write  this  as 


{E)  = 


WdE 


f-D 


+ 


LdE 


B  + 


WdE 


(6.10) 


(6.11) 


vr  Jw  VE2  -  W^'  ttJl  LVE2  -L2'  y/E‘^-W‘^'_  ' 

The  Hrst  and  third  terms  in  this  sum  have  the  same  integrand  and  so  can  be  combined: 

lr\-—  /~^  dU  ,  2L  dE 

~~  Jw  VE2  -  IU2'  Jl  VE2  -L2'  ■ 

These  integrals  are  inverse  hyperbolic  cosines  (ch~^):  the  general  identity  for  a  >  0  is 
f  dx 


(6.12) 


y/x‘^  —  a? 


=  ch 


-1 


=  In 


+ 


-  1 


With  this  identity,  (6.12)  becomes 


{E)  = 


vr 


VFch  ^:^  +  Lch 

W  L 


(6.13) 


(6.14) 


The  mean  extent  {E)  can  be  expressed  using  logarithms  by  way  of  (6.13).  Writing  the 
length-to-width  ratio  as  a  =  L/W ,  (6.14)  becomes 


w 


2 
vr  L 


In 


l-|-a2  -|-a)-|-aln(yl-|-l/a2  -|-l/a 


(6.15) 


Figure  23  plots  {E) /W  versus  a  from  (6.15),  and  shows  that  values  are  typically  around  2. 
That  these  values  are  reasonably  insensitive  to  the  ship’s  length  is  due  to  the  seeker’s 
narrow  beam,  as  discussed  in  Section  5.3. 
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Figure  23:  The  mean  extent  (E)  in  units  of  width  W  for  a  range  of  length- 
to-width  ratios  a,  from  (6.15) 


For  a  >  1,  the  term  in  brackets  in  (6.15)  is  approximately  ln(2a)  +  aln(l  +  1/a).  Now 
applying  the  identity  ln(l  +  x)  ~  x  for  |x|  ^  1  produces 

- — ^  =  —  (1  +  ln2  +  Ina)  «  (1.1  +  0.61na).  (6.16) 

W  TT 

So,  for  length-to- width  ratios  between  5  and  10,  {E)  ranges  from  2.1  IF  to  2.5  IF,  as  indeed 
is  shown  in  Figure  23.  But  measured  values  of  E  are  not  especially  peaked  about  (E) 
because  the  probability  density  p{E)  is  bi-modal. 


6.4  Cumulative  Probability  of  the  Extent,  C{E) 

The  cumulative  probability  of  measuring  any  extent  up  to  some  value  E  is 

rE 

C{E)=  /  p{E)dE.  (6.17) 

Jw 


In  principle  we  could  use  (6.9)  to  evaluate  this  integral  in  different  extent  regimes.  In 
practice  it’s  perhaps  easier  to  relate  p{E)  dE  to  2/7r  |d0|  as  we  did  when  calculating  p{E), 
which  amounts  to  a  change  of  variables  in  the  integration  in  (6.17). 


C{E)  for  W  ^  E  ^  L  (L-region  1) 

Integrating  from  IF  to  L  means  that  E  is  increasing  (i.e.  dE  >  0  so  |dE|  =  dE),  which 
corresponds  to  9  decreasing]  that  is,  d9  <  0.  Hence  |d0|  =  —d9,  and 


piE)  dE  =  p{E)  |dE|  =  2/7r  |d0|  =  -2/7rd0 
With  the  aid  of  (6.2),  the  integral  (6.17)  becomes 

r-E  o  /■  sin-1  (W/E) 

t/2 


c{E)=  r  —de=f 

Je=W  ^  Jtt 


—  de  =  — 

TT  7T 


sm 


-1 


IF  _  E 

"e  ~  2 


using  the  identity  sin  ^  x  +  cos  ^  x  =  tt  j 2  for  all  relevant  x. 


(6.18) 


—  cos 

TT 


IF 

"e  ’ 

(6.19) 
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Figure  24:  Cumulative  probability  density  C{E)  for  length-to-width  1.5, 
from  (6.24) 


C{E)  for  L  <  E  ^  D  (L-region  2) 

The  integration  (6.17)  now  requires  both  D-region  1  and  D-region  2.  Again  the  integration 
is  over  increasing  E,  so  dE  >  0,  meaning  |dE|  =  dE: 

p{E)  dE  =  p{E)  |dE|  =  2/7r  |d6»|£,_,egion  i  +  2/7r  |d(9|£,_region  2  •  (6-20) 

In  D-region  1,  increasing  E  means  increasing  9  (i.e.  d0  >  0  so  |d0|  =  d0).  In  D-region  2, 
increasing  E  means  decreasing  9  (i.e.  d0  <  0  so  |d0|  =  — d0).  Thus  (6.20)  becomes 


p{E)  dE  ‘Ij'K  d0£)_].ggjQn  \  ‘Ij'K  d0£)_].ggjQjj  2  • 
The  cumulative  probability  (6.17)  is  now  written 

C(E)  =  C(L)+  rp(E)dE  =  C(L)+  r  id0,)_,ggig,i-  r  - 
Jl  J  e=l  ^  Je=l  ^ 

which  becomes,  with  the  aid  of  (6.2), 

•cos~^{L/E)  2 


d9 


Z)-region  2  ? 


(6.21) 


(6.22) 


C{E)  =  C{L)  +  / 
Jo 


psin-^  (W/E)  2 

d9-  -d9 

^  J  sin-^  {W/L)  ^ 


2  /  _tW  ,  L 

=  —  cos  —  +  cos  — 

TT  V  E  E 


The  complete  cumulative  probability  is  then 


f  f  cos-1  w 

lU^  E 

C{E)  =  <^ 

if  (cos-1  f  +  cos-1  f) 

L  <  E  ^  D . 

(6.23) 


(6.24) 


This  is  plotted  in  Figure  24  for  a  ship  of  length  1.5  times  its  width.  (Again,  this  un¬ 
realistic  ratio  is  chosen  to  better  highlight  the  details  of  the  function.)  More  realistic 
cumulative  probabilities  are  plotted  in  Figure  25,  using  length-to-width  ratios  of  5  and  10 
respectively. 
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Length-to-width  ratio  5  Length-to-width  ratio  10 

Figure  25:  Cumulative  probabilities  from  (6.24)  for  realistic  length-to-width 
ratios.  The  rightmost  part  of  the  curve  becomes  less  significant  for  higher 
ratios 

As  a  concrete  example,  we  ask  for  C{{E)),  the  probability  of  measuring  any  extent 
less  than  the  mean,  for  the  case  of  a  =  L/W  =  5.  Here,  either  of  (6.15)  and  (6.16)  gives 
(E)  ~  2.10 VF.  Now  use  (6.24)  to  compute  C{2.10W),  referring  to  the  case  W ^  E  ^  L: 

C(2.10IF)  =  -  cos-1  ^  ~  0.68  .  ^g_25) 

TT  2.10 

That  is,  on  average  68%  of  the  extent  measurements  will  yield  values  between  W  and  2.10  W. 

Similarly  for  a  =  10  and  15  we  have  (with  the  above  a  =  5  result  included  for  com¬ 
pleteness) 


0  =  5 
o  =  10 
o  =  15 


(E)  ~  2.10 
(E)  ~  2.54 
(E)  ~  2.80 


and 

and 

and 


C{{E))  ~  0.68, 
C{{E))  ~  0.74, 
C{{E))  ~  0.77. 


(6.26) 


These  don’t  differ  markedly  from  the  corresponding  values  for  an  ellipse  with  narrow  beam 
in  (5.11). 


7  Concluding  Comments 

The  expressions  given  in  this  report  form  a  base  for  more  extensive  analyses  of  how  a 
seeker  might  interact  with  a  ship  based  solely  on  the  range  extent  measured  by  the  seeker. 

Any  analysis  must  use  a  tractable  model  of  the  ship  and  a  choice  of  beam  type  for  the 
seeker.  It’s  instructive  to  compare  the  statistics  for  the  elliptical  and  rectangular  models 
for  each  beam  type;  in  doing  so,  we  find  the  two  sets  of  statistics  to  be  similar  as  discussed 
in  the  next  paragraphs. 

For  the  wide  beam,  we  see  similar  values  for  extent  as  a  function  of  incoming  seeker 
angle  (Figures  2  and  8),  and  mean  extent  versus  length-to-width  ratio  (Figures  4  and  10). 
The  cumulative  probabilities  in  Figures  5  and  12  (left)  don’t  differ  greatly. 

For  the  narrow-beam  models,  the  extent  as  a  function  of  incoming  seeker  angle  is  just 
the  outline  of  the  ship:  Figures  14  and  20.  The  other  narrow-beam  results  cannot  be 
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compared  with  the  wide-beam  results,  but  within  the  two  narrow-beam  models  we  find 
that  the  two  probability  densities  versus  extent  are  similar  (Figures  15  and  22  left),  as  are 
the  mean  extents  versus  length-to- width  ratio  (Figures  16  and  23),  and  the  cumulative 
probabilities  (Figures  17  and  25  left). 

With  regard  to  computational  complexity  for  a  ship  model,  the  ellipse  is  generally 
more  straightforward  to  analyse  than  the  rectangle  because  the  ellipse  does  not  demand 
contributions  from  two  angle  bins  to  one  extent  bin,  unlike  the  rectangle;  this  is  a  very 
big  advantage  of  the  ellipse  model.  However,  the  integrations  required  for  the  ellipse  are 
generally  not  expressible  in  terms  of  elementary  functions,  unlike  those  for  the  rectangle. 
But  these  ellipse  integrals  are  easily  evaluated  numerically  so,  on  balance,  the  ellipse  model 
appears  to  be  the  better  choice. 

The  wide  beam  is  obviously  of  most  relevance  to  current  seekers,  but  this  could  well 
change  for  future  seekers.  Mean  ship  extents  for  the  wide  beam  are  shown  in  Figures  4 
and  10;  referring  to  (3.17)  and  (4.11)  with  the  comments  immediately  following  them,  we 
see  that  the  mean  extent  will  tend  to  have  a  value  near  the  ship’s  length,  especially  when 
the  length  is  much  larger  than  the  width.  In  contrast,  a  narrow  beam  produces  a  mean 
extent  closer  to  the  ship’s  width.  But  unfortunately,  the  situation  is  not  quite  as  simple 
as  this.  The  four  probability  densities,  ellipse  and  rectangle  wide  beam,  and  ellipse  and 
rectangle  narrow  beam  are  repeated  in  Figure  26  for  convenience,  and  show  that  p{E)  is 
mostly  bi-modal.  This  implies  that  in  any  one  measurement,  we  do  not  expect  to  measure 
the  mean  extent,  even  though  the  mean  extent  is  traditionally  called  the  extent’s  “expected 
value”.  Such  a  term  only  applies  when  p{E)  is  uni-modal.  Even  so.  Figure  26  indicates 
that  a  rough  rule  of  thumb  is  still 

-  wide  beams  tend  to  yield  extents  of  approximately  the  ship’s  length  (because  the 
probability  densities  have  more  support  at  the  right-hand  end  of  the  graphs) 

-  narrow  beams  tend  to  yield  extents  of  approximately  the  ship’s  width  (because  the 
probability  densities  have  more  support  at  the  left-hand  end  of  the  graphs) 

In  the  absence  of  any  other  information  except  beam  type  and  expected  ship  type, 
the  seeker  might  well  estimate  the  ship’s  size  by  assuming  that  a  value  “in  the  neighbour¬ 
hood  of”  the  mean  extent  has  been  returned  by  an  extent  measurement.  For  the  wide 
beam,  refer  to  Figures  4  and  10;  for  example  if  the  seeker  is  searching  for  a  ship  with  a 
length-to-width  ratio  of  10,  then  it  could  divide  its  extent  measurement  by  6  to  give  an 
estimate  of  the  ship’s  width.  (And  of  course,  the  length  would  then  follow:  L  =  lOVF.) 
Similarly  for  a  narrow  beam:  Figures  16  and  23  show  that  the  seeker  can  divide  its  extent 
measurement  by  about  2.4  to  give  an  estimate  of  the  ship’s  width.  But  this  idea  as¬ 
sumes  that  the  mean  extent  has  been  measured,  which  is  only  a  first  estimate  of  a  seeker’s 
strategy. 

The  seekers  in  this  report  approach  the  ship’s  centre  directly.  Having  a  seeker  approach 
off  to  one  side  would  not  give  very  different  results  for  a  wide  beam,  but  certainly  would 
yield  new  values  for  a  narrow  beam.  The  question  of  what  extents  a  narrow-beam  seeker 
would  measure  if  it  explored  different  look  angles  on  approach  would  extend  the  work  in 
this  report — but  does  become  very  model  dependent. 

In  the  final  analysis  the  seeker  sees  no  very  strongly  dominant  range  extents,  so  that  it’s 
probably  not  necessary  for  a  jammer  to  focus  heavily  on  constructing  any  one  particular 
range  extent  when  assembling  a  false  target  presentable  to  a  seeker.  Of  course,  a  caveat 
is  that  measuring  a  ship’s  radial  extent  is  just  one  strategy  that  a  seeker  might  employ. 
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Ellipse,  wide  beam;  L/W  =  5  (Figure  3)  Rectangle,  wide  beam;  L/W  =  10  (Figure  9) 


Ellipse,  narrow  beam;  L/W  =  5  (Figure  15)  Rectangle,  narrow  beam;  L/W  =  10  (Figure  22) 
Figure  26:  All  probability  densities  repeated  here  for  convenience 


The  results  of  this  report  allow  a  selection  of  more  representative  range  extent  values 
for  ongoing  work  into  false-target  generation  for  countering  incoming  missiles.  The  main 
result  here  is  that  they  indicate  we  have  no  strong  obligation  to  manufacture  a  waveform 
that  mimics  the  full  length  of  a  ship.  This  allows  more  economical  use  to  be  made  of 
available  false-target  generation  resources. 

I  wish  to  thank  Dr  Roland  Keir  of  CEWD  for  his  thorough  suggestions  regarding  the 
content  and  direction  of  this  report. 
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